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C$ . Abstract 

We analyze the finite-size corrections to the energy and specific heat of the 
critical two-dimensional spin- 1/2 Ising model on a torus. We extend the analysis 
of Ferdinand and Fisher to compute the correction of order L~ 3 to the energy 
and the corrections of order L~ 2 and L~ 3 to the specific heat. We also obtain 
general results on the form of the finite-size corrections to these quantities: only 
integer powers of L _1 occur, unmodified by logarithms (except of course for the 
leading logL term in the specific heat); and the energy expansion contains only 
odd powers of L . In the specific-heat expansion any power of L _1 can appear, 



but the coefficients of the odd powers are proportional to the corresponding 



coefficients of the energy expansion. 
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1 Introduction 

It is well-known that phase transitions in statistical-mechanical systems can occur 
only in the infinite- volume limit. In any finite system, all thermodynamic quantities 
(such as the magnetic susceptibility and the specific heat) are analytic functions of 
all parameters (such as the temperature and the magnetic field); but near a critical 
point they display peaks whose height increases and whose width decreases as the 
volume N = L d grows, yielding the critical singularities in the limit L — > oo. For 
bulk experimental systems (containing N ~ 10 23 particles) the finite-size rounding 
of the phase transition is usually beyond the experimental resolution; but in Monte 
Carlo simulations (N ^ 10 6 -10 7 ) it is visible and is often the dominant effect. 

Finite-size scaling theory [|], 0, Q, || provides a systematic framework for under- 
standing finite-size effects near a critical point. The idea is simple: the only two 
relevant length scales are the system linear size L and the correlation length ^ of 
the bulk system at the same parameters, so everything is controlled by the single ratio 
£oo/L.[] If L 3> £oo, then finite-size effects are negligible; for L ~ £00, thermodynamic 
singularities are rounded and obey a scaling Ansatz O ~ L po Fo(^ OQ /L) where po is 
a critical exponent and F@ is a scaling function. Finite-size scaling is the basis of the 
powerful phenomenological renormalization group method (see Ref. || for a review); 
and it is an efficient tool for extrapolating finite-size data coming from Monte Carlo 
simulations so as to obtain accurate results on critical exponents, universal amplitude 
ratios and subleading exponents [|5], [| ^, |8|, and references therein] .0 In particular, 
in systems with multiplicative and/or additive logarithmic corrections (as the two- 
dimensional 4-state Potts model [TO]), a good understanding of finite-size effects is 
crucial for obtaining reliable estimates of the physically interesting quantities. 

In finite-size-scaling theory for systems with periodic boundary conditions, three 
simplifying assumptions have frequently been made: 

(a) The regular part of the free energy, / reg , is independent of the lattice size L || 
(except possibly for terms that are exponentially small in L). 

(b) The scaling fields associated to the temperature and magnetic field (i.e., g t and 



ghi respectively) are independent of L ]TT|. 



(c) The scaling field gi associated to the lattice size equals L 1 exactly, with no 
corrections L~ 2 , L™ 3 , ... f§]. 

Moreover, in the nearest-neighbor spin-1/2 2D Ising model, it has further been as- 
sumed that 



(d) There are no irrelevant operators [12], [T3 



1 This is true only for systems below the upper critical dimension d c . For Ising models with 
short-range interaction, d c = 4. 

2 Finite-size scaling has also been successfully applied to data coming from transfer-matrix com- 
putations H. 



Unfortunately, the combination of these four assumptions implies that the asymp- 
totic expansions for the energy and specific heat for the Ising model at criticality 
terminate at order 1/L (see Ref. ||14|| ). However, the numerical results presented in 



||14] , [7j (as well as the analytic results presented in this paper) show this to be false. 
The problem, therefore, is to determine which one(s) of these assumptions are in- 
valid, and why. Assumption (c) is extremely plausible from renormalization-group 
considerations, at least for periodic boundary conditions; and assumption (d) has 
been confirmed numerically through order (T — T c ) 3 at least as regards the bulk be- 
havior of the susceptibility |K|. However, both numerical [|l|, [16! anc ^ theoretical 117 1 



evidence has recently emerged suggesting that irrelevant operators do contribute to 
the susceptibility at order (T — T c ) 4 . 



In a classic paper, Ferdinand and Fisher [18|] considered the energy and the specific 
heat of the two-dimensional Ising model on a torus of length L and aspect ratio p, 
and obtained the first two (resp. three) terms of the large- L asymptotic expansion of 
the energy (resp. specific heat) at fixed x = L(T — T c ) [this is the finite-size-scaling 
regime] and fixed p.f\ In particular, at criticality (T = T c ) they computed the finite- 
size corrections to both quantities to order L" 1 . In 1999, Hu et al. |22|] published 
(without details) the correction of order L~ 3 to the energy and showed that the L~ 2 
correction is absent .0 In this paper we compute explicitly the following finite-size 
corrections: 

• The correction of order L~ z to the energy. 

• The corrections of order L~ 2 and L~ 3 to the specific heat. 

Furthermore, we also find new insights into the general analytic structure of the 
finite-size corrections to this model. We show that in the critical two-dimensional 
Ising model: 

• The finite-size corrections to the energy and specific heat are always integer 
powers of L~ l , unmodified by logarithms (except of course for the leading logL 
term in the specific heat). 

• In the finite-size expansion of the energy, only odd integer powers of L _1 occur. 

• In the finite-size expansion of the specific heat, any integer powers of L~ x can 
occur. However, the coefficients of the odd powers of L _1 in this expansion are 
proportional to the corresponding coefficients in the energy expansion. 



3 Ferdinand and Fisher JTa] also obtained the position of the maximum of the specific heat Xmax (p) 
which depends on the torus aspect ratio p. Furthermore, Kleban and Akinci J19, E0[ showed that an 
excellent approximation can be obtained by keeping only the two largest eigenvalues of the transfer 
matrix, and they interpreted their results in terms of domain-wall energies. This approximation 
is already good at p = 1 and becomes exponentially better with increasing p. This method could 
surely simplify the computations presented in this paper; but here we are interested in the exact 
values of the finite-size-scaling corrections. These results are used as theoretical inputs in H, Elf. 
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The published version of [E2| contains several misprints in the crucial formula (3.23). Version 



2 of this paper in the Los Alamos preprint archive cond-mat contains the correct formula. 



These results can be compared to the general renormalization-group expression for 



the finite-size corrections to the energy and the specific heat |L4[], in which arbitrary 
powers of L~ l and terms of the type L _1+p logL (where p is some real number) 
can occur. The implications of these results for understanding which one(s) of the 
assumptions (a)-(d) are invalid will be analyzed elsewhere [pTf . 

The plan of this paper is as follows: In Section |2] we present our definitions and 
notation (generally following Ref. |H|). In Sections |3] and f| we present the computa- 
tion of the next terms in the asymptotic expansions for the energy and specific heat, 
respectively. Finally, in Section |5] we present our arguments about the type of finite- 
size-scaling corrections that can occur in these two expansions. We have summarized 
in Appendix |A| the basic definitions and properties of the ^-functions that will be 
needed in this paper. In Appendix [TJ we recall the Euler-MacLaurin formula. 

Remark. After the completion of this work, we learned that similar results have 
been independently obtained by Izmailian and Hu 



2 Basic definitions 

Let us consider an Ising model on a torus of size mxnat zero magnetic field. 
The Hamiltonian is given byQ 

H = -K^2a i a j (2.1) 

<<J) 

The partition function can be written as 

i 4 

Zmn = £V W = -(2 Sillh 2K ) mH/2 J2 Zi ( K > n ' m ) ' ( 2 ' 2 ) 

M i=l 

where the partial partition function Z\ is given by 

71-1 

Z 1 (^n,m) = n 2cosh (^T ±1 ) (2 ' 3) 

r=0 

and the rest are defined analogously using the first three columns of the following 
table [the last two columns will be needed afterwards in (2.9)/ (|4.18 )] . 



Zx 

Z 2 

z ?t 
z± 



2r + 1 cosh tanh sech 

2r + 1 sinh coth i csch , . 

2r cosh tanh sech 

2r sinh coth i csch 



5 In this paper we are following basically the notation used by Ferdinand and Fisher in |Q , with 
a few minor modifications. 



The quantities 7/ = ji(K, n) are defined by 



In 

cosh 7/ = ci = cosh 2K coth 2K — cos I — 

n 



In particular, we have 



7o 

7; 



2K + \og{tajih.K) 



log( Q+ \/cf -1 



1^0 



(2.5) 

(2.6a) 
(2.6b) 



The quantities 7; (2.6) satisfy 7; = ^n-i, and 7/ is a monotonically increasing function 
of I for < / < n. 

The internal energy density E and the specific heat Ch are given by 



E(K,m,n) 



coth2if 



mn 





4 



(2.7) 



C H (K,m,n) 



-2csch. 2 2K + 



E^' 


i=l 


2" 


4 
E^ 

_i=l 


4 

\§ z v 





(2.8) 



where the primes denote derivatives with respect to the coupling constant K. In 
computing observables (2"77)/ Q2~ffi ), the following formulae, derived from fl2.3| ), will be 
useful: 



n-l 



Yi = yL^ +1 tanh(— 



Z x 



r=0 
n-l 



(2.9a) 



— 2^7 2 r+itanh^ — 



r=0 



+ 



n-l 



171 V^ 'I V U /^ m 72r+l 

y2^ 7 2r+itanh^ — - — 



r=0 



in 



r=0 



(2.9b) 



The analogous ratios for i = 2, 3, 4 can be obtained from (|2.9aj) /(2.9b) by using (|2.4j ) 



[Note that the third column of ( |2.4j ) does not play any role here]. The factor i in the 
entry i csch of ( |2.4|) changes the sign of the last term of (2.9b) for Z'% and Z". 

The critical point of the Ising model corresponds to the self-dual point sinh 2K C = 
1. That is 

K c = 1 log(l + y/2) (2.10) 

In this paper we are concerned with the finite-size-scaling corrections to the energy 
and specific heat of the critical Ising model. We will express all our results in terms 



of the length n and the aspect ratio of the torus p:[] 



m 

P= — 

n 



(2-11) 



Indeed, all our results are invariant under the transformation n <-> m. Here we shall 
show that the energy and specific heat at criticality have asymptotic expansions of 
the form 



- E(K c ,m,n) = -E c (n, p) = E +) j 



E k {p) 



fc=i 



rr 



C H (K c ,n,m) = C H , c {n, p) = C 00 \ogn + C (p) + ^ 



C k (p) 



fe=i 



tr 



(2.12) 
(2.13) 



The coefficients E and Coo can be obtained from Onsager's solution [25fl; Ei, Cq, and 
Ci were computed by Ferdinand and Fisher |18| ; and finally, the fact that E 2 = and 
the expression for E 3 were given (without details) in Ref . [Q . Here we shall compute 
explicitly the terms E 3 , C 2 and C3, and shall show that E 2 = E4 = E G = . . . = 0. 

Let us now see how 7; and its derivatives behave close to the critical point ( |2.1(J| ). 
To do so, we introduce the finite-size-scaling parameter r as in Ref. [JT8|, Eq. (2.12)][| 



T\2 
71/ 



- lsmh2K+ — 

2 V sinh2ir 



(2.14) 



Thus, t = corresponds to the critical point K = K c , and r^O fixed corresponds 
to the finite-size-scaling regime n — > 00, K -^ K c with n(fr — K c ) fixed. Hereafter, 
we will consider the behavior of all quantities as a function of r in the limit r — > 0. 
The value of 70 at r = is zero; its behavior close to the critical point is given by 

70 (r,n) = - 



nJ 


w\ 



The derivatives of 70 with respect to K are non- vanishing at criticality: 

ofTo =4 

dK 

d 2 7o 



To (0,n) = 
7o(0,n) = 



T=T C 



dK 2 



-4V2 



(2.15) 

(2.16a) 
(2.16b) 



T=T C 



(Note that prime continues to denote d/dK, not d/dr. However, the final result will 
be expressed in terms of r [in the limit r —>■ 0], hence the notation 7 (0,n)). For a 
generic / 7^ the critical value of 7/ is given by 



7 *(0,n)=21og 




+ Sln (£) 



(2.17) 



where 



6 In conformal-ficld-theory language, the modular parameter of the t orus i s t = i p 
this t has nothing to do with the temperature- like parameter defined in (|2.14j) . 

7 The parameter r plays the same role as the usual finite-size-scaling parameter x = L(T — T c ). 
Indeed, to leading order in K — K c we have r = —2n(K — K c ). 
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while its derivatives with respect to K are given by 

c' 



7j 



ri 



Vcf^i 



ci{c[) 



,/\2 



V^i (<* - !) 3/2 

where the quantity c[ vanishes at criticality as 



(2.18a) 
(2.18b) 



c^r^n) =c[r,n) 



®+° 


[01 



and the quantity d! gives a non-zero value 



c"(0, n) = c"(0, n) = 16 



(2.19) 



(2.20) 



We can write the partial partition functions Z{ [cf. (|2.3| )/ (|2.4j )1 in the following 
form: 



n-l 



2r+l 



r=0 



Z x (r, n, p) = Pi (r, n, p) exp I — ^ 7 2r+ i 

(n-l 
y5Z 72r 

Z 3 (r, n, p) = P 3 (r, n, p) exp f ^ ^ 7 2r j [l + e""^] 

Z 4 (r, n, p) = F 4 (r, n, p) exp f ^ ]T 72r ) [l - e' m ^] 
where the quantities Pj(r, n, p) are given by 

n-l 

log Pi (r,n,p) = ^log(l + e" 

r=0 
n-l 

log P 2 (r, n,p) = Y^ 1 °g( 1 ~ e ~" 

r=0 

n— 1 

logP 3 (r,n,p) = ^log(l + e" 

r=l 

n-l 

logP 4 (r,n,p) = ]Tlog(l-e- m ^) 



-ni72r+l 



-m72r+l " 



-ni72r N 



r=l 



The functions ( |2.22| ) give non-vanishing constants in the limit r — > [|P8 

]ogP x (0,n,p) = ^ + 0(n- 2 ) 



(2.21a) 
(2.21b) 
(2.21c) 
(2.21d) 



(2.22a) 
(2.22b) 
(2.22c) 
(2.22d) 



(2.23a) 



logP 2 (0,n,p) = 6 -± + 0(n- 2 ) (2.23b) 

logP 3 (0,n,p) = \ e -l e ^ 4 + 0(n- 2 ) (2.23c) 

logP 4 (0,n,p) = 9 2 + O(n- 2 ) (2.23d) 

where the functions 9i with i = 2,3,4 are the usual ^-functions (see Appendix |A|), 
and 9q is defined in ( |A.2| ). 

Finally, we introduce the ratios 

Zi(r,n,p) 



Ri(T,n,p) 



Zi{r,n,p) 



Thus, from ( |2.21| ) we get 

Ri(r,n,p) 
R 2 (T,n,p) 



P2(r,n,p) 
Pi(r,n,p) 



R 3 {T,n,p) = 2 cosh — - )P (T,n,p) — 

V 2 / P^r.n.p) 

R 4 {T,n,p) = 2smh — - P (r,n,p) — ^ 

V 2 / P x (T,n,p) 



where Po(T,n,p) is defined as 



log P (r,n,p) 



r/i 



n-l 



n-1 



X)T2t~X)^ 



r+1 



r=l 



r=0 



The sum of the four ratios is denoted by R 



R(r,n,p) = ^i?i(r,n,p) 



4=1 



The ratios P 2 and P3 have a non-vanishing value at the critical point ||T 

R 2 (0,n,p) = °^ + 0{n- 2 ) 

R 3 (0,n,p) = ^ + 0(n- 2 ) 
#3 

while R4 vanishes at K = K c : 

R 4 (T,n,p) = -sinh(rp) [0 2 4 + 0(rr 2 )] + C(r 2 
The sum of the four ratios at criticality is a non-zero constant 

P(0,n,p) = + 0{n ) 

The function P has also a non-vanishing limit at criticality: 

P (0,n,p) = e- 7Tp/4 [l + O(n 



(2.24) 

(2.25a) 
(2.25b) 

(2.25c) 

(2.25d) 

(2.26) 

(2.27) 

(2.28a) 
(2.28b) 

(2.29) 

(2.30) 

(2.31) 



3 Finite-size-scaling corrections to the internal en- 
ergy 

The internal energy at the critical point E c is equal to 

1 4 Z' 
-E c (n,p) = -E(K c ,n,p) = V2 + lim—^J2f R * t 3 ' 1 ) 

The terms Z^/Zi with % = 1, 2 vanish trivially as all the 7 2r+ i vanish. The term Z'^jZ^ 
does not vanish due to the contribution of 7q; but its total contribution is also zero as 
it is multiplied by tanh(m7o/2), which vanishes at criticality. The only non- vanishing 
contribution comes from % = 4: 



So we obtain the formula 



_i = ^ C othf^) ~ -2mcoth(pr) as r^O (3.2) 

Za 2 \ 2 / 



E c (n, p) = V2-- lim ^'^ coth(pr) 
n t-*o R(0,n, p) 



(3.3) 



where _R(0, n, p) is given by ( |2.30[ ). Note that, by fl2.29| ), Ra.(t, n, p) ~ r as r — > 0, so 



.R 4 (t, n, p) coth(pr) gives rise to a non-zero result in this limit. 

The goal of this section is to extend the Ferdinand- Fisher asymptotic expansion 
[|T8| to order n~ 4 . Let us first consider the quantity logP 4 at criticality: 

n— 1 co . L"2"J 

lo g P 4 (0,n,p) = ^log(l- e - m ^) =-2j2~J2e- m ^ (3.4) 

r=l p=l r=l 

where [^J is the largest integer < x. The sum over r can be split into two parts: 

LfJ »(«)-i LfJ 

£=£ + £ ( 3 - 5 ) 

r=l r=l r=s(n) 

where s(n) will be chosen afterwards. We can drop the second sum in (|3.5|) , as it 
gives a contribution of order ~ n\ log[l — exp(— m72 S )]| ~ nexp[— 2irps(n)]. Instead 
of the choice s(n) = (3/2irp) logn made by Ferdinand and Fisher, we shall use the 
choice 

s(n) — - — logn (3.6) 

27rp 

with M an arbitrary positive integer, to ensure that the total contribution of the 
second sum in ( |3.5|) is as small as we want (namely, ~ n - ^ -1 )). 
We can use the following expression for 7 2r at criticality fl"8|| : 



-7 2r (0,n) = log 



s " l, ir) + V 1 + sm2 



(3.7) 



to obtain an asymptotic series of m^r in terms of n . 

m l2r (0, n) = 2npr - ^4 + ^4 + °( n ~ & ) ( 3 - 8 ) 

3 n z on 

We should recall that the choice of s(n) ([3.6|) for any integer M guarantees that the 
ratio r/n is a small quantity for < r < s(n) (as r/n < s(n)/n <C 1 if n is large 
enough) . 

Remark: It is interesting to note that only even powers of n~ x occur in the expansion 



Plugging the expansion ( |3.8| ) in Eq. ( |3.4| ) we obtain 

-2vrrpp *" P 

/ ° 



oo 1 s(n)— 1 

logP 4 (0,n,p) = -2^- ^ e 



3 oo s(n)-l .j 

rl E ' 

p=l r=l p=l r=l 



E E ^ e " 27rrpp + °^~ 4 ) ( 3 - 9 ) 



We can extend the sums X/r=i ^o S^=i m (P-^D by introducing an error of order 
n -M_ fj^g secon( j term of the r.h.s. of ( |3.9| ) can be expressed in terms of 



E 

r=l 



r 3g-27rrpp 



sinh (7rpp) 2 sinh (irpp) _ 



Finally, we can write P&{n, p) in the following form 



P 4 (0,n,p)=9 2 



i-l*jfW 



which improves ( |2.23d| ) . The function pi (p) is given by 

r 3, °° r 



P1(P) 



2vr J p 



E 

m=l 



sinh (rmrp) 2 sinh (rmrp) 



Using similar methods one can obtain the improved version of ( 2.23 ): 



Pi(0,n,p) 

A(0,n,p) 
A(0, n,p) 



E: 
Ei 

\e_ 

26 



1 



Pi 



1 + _L(|_ P2 ) + O( o 

1 pi 



2 g irp/4 



where p 2 , pi, and p 2 are defined by 

27T d p 



Pi(p) 

P2(p) 
P2(P) 



E(-i) m+1 



m=l 

3 °° 



1 + _ ¥ + 0(0 



1 



EE V^ 

9/L 2_^ 



m=l 



24 

Q oo 

irp 

~2A 



sinh (mirp) 2 sinh (mirp) 
3 1 



E(-d 



sinh (mirp/2) 2 sinh (mirp/2) 
1 3 



m+l 



m=l 



1 



sinh (miip/2) 2 sinh (mirp/ 2) 
10 



(3.10) 

(3.11) 

(3.12) 

(3.13a) 
(3.13b) 
(3.13c) 

(3.14a) 
(3.14b) 
(3.14c) 



Finally, we have to improve the expression of log Pq ( |2.31|) . Let us first consider 
the sum 



rn 
~2 



n-l 



n-l 



5^72r(0,n) = m^log 



r=l 



r=0 



fTTX 

sin — 
V n 



!fsin 2 (^) 
n / 



(3.15) 



We can apply the Euler-MacLaurin formula ( |B.4|) to the function f(p) = log(sinp + 



a/1 + sin 2 p) with L = 2n and a = 1/2. The result is^ 

2mn Tip 7T 3 p 1 



n-l 



m 

~2 



$^72,(0, 



n 



-G 



r=l 



7T 



6 180 n 2 



+ C(n" 



(3.16) 



where G ~ 0.915965594177219 is Catalan's constant. Using similar methods we 
obtain the other sum appearing in ( |2.31| ): f\ 



m 

~2 



n-l 



£ 7a . +1 (0,n) = — G 



2mn i\p 7ir 3 p 1 



r=0 



12 1440 n 2 



C(n" 4 ) 



(3.17) 



Putting ( p,16| ) and (|3.17|) together we obtain the improved version of ( |2.31|) 



P (0,n,p)=e-^ 
where pz(p) is defined as 



Ps(p) 



11 r 



+ 0(n 



7T" 



Ps(p) = ^P 



(3.18) 



(3.19) 



The improved expressions for the ratios Ri(r = 0,n,p) are easily obtained from 

(EH/dH) 



R 2 (Q,n,p) = -£- 
R 3 (0,n,p) = -£- 

#3 



1 l I Pi , ~ Pi \ , , -i, 



rr 



1 + -a (Pi - P2 - Ps) + 0{n 



(3.20a) 
(3.20b) 



_R 4 (0,n, p) = — sinh(pr)#2#4 



+ 0{n 



i l fPi , ~ Pi , 

1-^ y+P2- y +P3 



(3.20c) 



Plugging the formulae ( |3.20| ) in the expression for the critical energy density (|3.1| ) 
we get 

-E c (n,») = ^ + ^l + m + o(l i ) (3.21) 

8 It is easy to verify that f^\p) is integrable over [0, it]. This means that the next term in the 
expansion (3.16) is of order 0(n~ 4 ). 

9 The leading terms of Eqs. ( 3.16|) /( 3.17| ) were obtained by Ferdinand |2q| . 
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where E\(p) |1| and E 3 (p) are given by the expressions 



E ^ = e^iTT, (3 ' 22) 

E3{p) = ~ {e 2 + lT+\r { Mp) i 04 + 2) - p2{p)6 ' 

+ Pi(p) (^y^) +Up)0z+Vz(p)(0; + 6^ (3.23) 
The numerical values of the function E 3 (p) are given in Table [1] 



P E 3 (p) 



1 -0.206683145336864 

2 -0.184202899115749 

3 -0.153247694215529 

4 -0.102599506933675 

5 -0.061201301359728 

6 -0.034200082347112 

7 -0.018369506074164 

8 -0.009614465215356 

9 -0.004941568941314 
10 -0.002505707497764 
15 -0.000074110828658 
20 -0.000001946957522 
00 



Table 1: Values of the coefficient E 3 (p) from ( |3.23|) for several values of the torus 
aspect ratio p. 



Remarks. 1. After the completion of this work, Prof. Izmailian informed us that the 
correct expression for the coefficient E 3 (p) had been published in the revised version 
of Ref. |Z5| (which can be found in the Los Alamos preprint archive cond-mat). Their 
expression is surprisingly simple 

^> = -^OT^^] (3.24, 

Indeed, the numerical value of ( p.24[ ) coincides with our result ( |3.23| ). It would be 
interesting to find the analytic identities proving the equivalence of ( |3.23|) and Q3.24J ). 
2. Let us check that ( |3.23| ) has the correct behavior under m <-> n (p <-> 1/p). 
Indeed, the (trivial) fact that E c (n,m) = E c (m,n) implies that we should have 

Bk{p) = ^^ (3.25a) 

P 

E 3 (p) = ^M (3.25b) 

P 3 

12 



The first equation (|3.25a|) can be easily proved by using Jacobi's imaginary transfor- 
mation of the ^-functions ( |A.4|) . Using these transformation one can easily show that 
the second equation (3.25b) holds for ( |3.24| ). We have also verified numerically that 



(3.25b) holds for our result ( |3.23|) to high accuracy using mathematica. 

3. There is another simple way to test our results: We can first compute the exact 
value of the critical energy density E c (n, p) for several values of n and a fixed value of 



p by using fl2.7| )/( |2T3"D /( |2.4D . Then, we subtract the first two terms of the expansion 



( p.21| ) and fit the resulting function to the Ansatz B 3 n~ 3 + B 5 n~ 5 + B 7 n~ 7 . In Table |2] 
we show the numerical results for such fits with p = 1,2, and 3. For p = 1 we have 
used in the fits 1309 different values between n = 16 and n = 4096. For p = 2, 3 
we have used 10 different values corresponding to n = 2 P , p = 2, . . . , 11 (that is why 
our estimates are more accurate for p = 1 than for p = 2,3). We find an excellent 
agreement among the numerical estimates for B 3 and the exact values of E 3 quoted 
in Table [I]. The value of B 5 for p = 1 also agrees well with the value ~ —0.7301823 



obtained by Izmailian [£7| using analytic means 



B 3 -0.2066831453369 -0.1842028991157 -0.1532476942155 
B 5 -0.73018231235 -0.416996817 -0.316738073 

B 7 -4.9362 -3.405 -2.693 



Table 2: Fits of the function -E c {n,p) - \/2 - E x {p)/n [cf. (|3~2^) 1 to the Ansatz 
B 3 n~ 3 + B 5 n~ 5 + B 7 n~ 7 for several values of the torus aspect ratio p. 



4. In the limit p —>■ oo of an infinitely long torus (i.e. a cylinder) we have 

E t (oo) = E 3 (oo) = (3.26) 

as lim^oo p6 2 = 0, and lim^oo 6> 3 = lim^oo 6 4 = 1 [cf. ( |A.1|) ]. 



4 Finite-size-scaling corrections to the specific heat 

The goal of this section is to extend the asymptotic series of Ferdinand and Fisher 
18| for the specific heat through order n~ 3 . Let us start with the definition (|2.8| ) 
and see which terms contribute to the critical value of Ch- The first term in ( |2.8j ) 
is just a constant (= —2); while the third term is quite similar to the one already 
obtained for the energy density (|3T3| ) [= — ApR~ 2 R\ coth 2 (rp)]. The most involved 
term is the second one. Using the analysis of Ferdinand and Fisher, we can get the 
final expression for the critical specific heat Cn,c(n,p) = Cn(K c ,n, p): 

n i \ n , /n 4i? 3 (0,n,p) R 3 {0,n,p) 

R(0,n,p) R(0,n,p) 

13 



R(Q,n,p)j^ 



) Ri(0, n, p)Qi ti H hm — — coth rp 

^—? n t^o R(u, n, p) 

(4.1) 



-4p lim I — ; — - — ^— — coth rp 



where the Qi } ± and Qi,i are those defined in Ref. [jn| evaluated at r = 0: 

1 _ tanh (Hn|±i) 



Q 



i.i 



g 



1,2 



Ql,3 
Ql,4 

Qi,- 



n,p) 

n,p) 

n,p) 
n,p) 

n,p) 

n,p) 



1 71—1 



r=o sin 

n-l 



(r+l/2)7r 



1 + sin' 



2 f (r+l/2)ff 



1/2 



E 



l - coth (nn|±i) 



»-=o sin i '—>— 

' n 



1 + sin' 



2 ((r+l/Zyn 



1/2 



n-l 

-E 



1 - tanh (^) 



sin(^)[l + sin 2 (^)] 1/2 
1 - coth (23&) 



r=l 
n-l 

S^n(=)[lW( ? )]" 2 

n-l 

-E 



n ^o ■ ( (r+l/2)7T 

r=o sm L ^- 

V n 

n-l 



1 + sin 



2 f (r+l/2)7r 
n 



1/2 



n ^ 



sin(^)[l + sin 2 (^)] 



1/2 



(4.2a) 

(4.2b) 

(4.2c) 
(4.2d) 
(4.2e) 

(4.2f) 



The terms with the factors Ri/R (i = 3, 4) can be obtained easily using the results 
of Section |3|. Let us first consider the quantity Qi >+ (n, p) ( |4.2f| ). The first step consists 
in expanding the factor [1 + sin 2 (r7r/ri)]~ 1 / 2 in Eq. ( |4.21| ) in power series of sin(r7r/ri): 



n-l 



n— 1 oo 



Ql.H 



n^sin(^) n^-f£-f 

r=\ \ n J r=\ k=\ 



T)^iT 



)« 



(2) 



= q\x+q\: 



(4.3) 
(4.4) 



The computation of Q\ + is done by applying the Euler-MacLaurin formula ( |B.4|) to 
the function f(p) = sin 2fc_1 (p) with L = 2n and a = 1/2. The result is 

log 2 n 1 



QT+(n,p) 



n + 12 n 2 



+ 0(n 



(4.5) 



The computation of the divergent part Q 1+ is & little more involved. The idea if to 
apply the Euler-MacLaurin formula ( P-4j) to the function f(p) = sin _1 (p) — 1/p + 
l/(j) — 7i") with L = 2n and a = 1/2: 



1 n— 1 

i E 

r=0 



, /m\ n 
sin I — i 



n 



\ n 7 tit 7r(r — n) 



-i°«- + A(A-s) ■«""" 
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ra-1 



)(D 



r=l 

Using the well-known asymptotic expansion |2!| [See also Q5.5D 1 



(4.6) 



-4^ 



r=l 



(4.7) 



(where 7s ~ 0.5772156649 is the Euler constant) we finally get 



TV 



. 2 7T 2 

log n + 7s + log - 



vr 72n 2 



C(n- 4 ) 



(4.8) 



Putting together fl4.5|) /(|Pj) we arrive at the final result 



Qi,+ (n,p) = - 



n 



2 l/2 



7T" 



log n + 7s + log h — -=■ + 0(n 4 ) 

7r 36rr 



(4.9) 



Using similar methods we obtainF] 

2 



Qi,-(n,p) 



71 



2 5/2 n 2 

log n + 7^ + log — -^^ + 0(n 



IX 



72n 2 



(4.10) 



The last part consists in evaluation the Qi,j (with i — 1, 2, 3) in (|4.2j ). For brevity 
we will do explicitly the simplest case Qi^(n,p) ( |4.2d|) . The first step is to expand 
the term 1 — coth(m7 2r /2) as a power series in exp(— m^r) 



- 

L -, 1 



2 J OO 






-mp72r 



=^=^in(^)Jl + sin 2 (^) 



(4.11) 



Then we split the sum over r as in (|3.5|) . The term including the sum Ylr=s(n) gi ves 



a total contribution of order n ( M ^ with the choice (|3.6|) for s(n). The second step 
is to plug into (|4.11|) the expansion (p.8[) for m7 2r 



00 s— 1 



QlA 



n 2-^2-^ 



-2irrpp 



i-^in^JH-sin 2 ^) 



_8vr 3 p 
3n 3 

Vl,4 ~r Wl,4 



-2nrpp 



00 s— 1 

P =/ r =i r sin(^) x /l-:-,ur(--) 



+ C»(n- 4 ) (4.12a) 



(4.12b) 
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The leading term of this equation was obtained by Onsager [E5j 
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The computation of Q\ \ is quite easy: we expand the factors sin(r7r/n) in powers 
of rir/n: 

2 °° s(n)-l 
p=l r=l 



Then we can extend the sum ^21=1 to S^Li a ^ an error of order n ( M + 2 ) log 
Using the fact that 



n. 



00 p — 2-KTp 1 

E ve -2nrpp = f = 1 

_ x (i_ e -2^p) 2 4sinh 2 (7rrp) 

we have that 

# > = _!^y;_^_ ( 4.i5) 

3n 2 ^ S i n h 2 (7rrp) 
The computation of Qi i follows the same steps: 



00 00 ^ . 00 00 






7T * — ' z — ' r 3n 2 

p=l r=l p=l r=l 



In this case the error introduced by extending the sum YHr=l to Y^T=i ls °f or der 



s(n)-l 



/? 



-ill 



. Using (|4.14|) we find that 



Q<°> = - V -[I - coth(Trrp)] -p-Yr[l- coth(Trrp)] + 0( n - 4 ) (4.17) 

7r z — ' r 3n z z — ' 

Thus, we write the final result as 

{2) (o) 
Qu(n,p) = Qfl(p) + ^±4p + 0(n" 4 ) (4.18a) 

n 

2 °° 1 
QS5W = -$>[l-coth(7rrp)] (4.18b) 

r=l 
27T f °° °° 1 

gj 2 J(p) = <J ^V[l -coth(vrrp)] + vrp^r 2 [icsch(7rrp)] 2 }> (4.18c) 



r=l r=l 



The other three quantities Qi,i (|4.2a| -c) can be computed in a similar way. The result 
can be written as ( |4.18[ ) using the translations given by ( |2.4j) [In this case, the third 



column of ( |2.4| ) does not play any role]. For i = 1,2 we should make two slight 
modifications: (a) The factor r in (|4.18|) should be replaced by r + 1/2, and (b) the 
sums over r in Qi,i an d Qi,2 start at r = rather than at r = 1 (as in Q 13 and Qm)- 
Putting all the pieces together we arrive at the final result 

C HA n, p) = I log „ + C„(„) + ^lM + ^ + ^ + O (A) (4-19) 
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where the coefficients Ci(p) are given by 

8 (, 25/2 
c o(P) = - log + 1e - 7 

71 \ 71 4 



Ci{p) 

c 2 (p) 



-2^ 



$2 + $3 + $4 
#2 #3 #4 



71 *-^ 



P 



a2a2n2 

V 2 V 3 V 4 



u=2 



@2 + $3 + #4 



-4p 



#2 + #3 + $4 
620364, 



•V2Ei{p) 



(4.20) 
(4.21) 



^2 + ^3 + ^4 E3(P) -^- 4A3(/?) ^ + L^ 



log 16 



7[ 



71 



02 



lfi £Mog0„ 

G(p) 



C 3 (p) 



#2 + $3 + $4 
•V2E 3 (p) 



where 



G(p) 

4i(p) 

AM 



)( 2 )/ 



>(2)/ 



)(2) 



(0) 



)(0) 



Qi> + Q$0 4 + Q\%6 2 - Q^A 2 (p)9 4 - QZiA 3 (p)6 [ 

P2(P)+P2{P) - 2 [P1(P) +P1(P)1 
P3(P)+P2(P) -Pl(p) 



(4.22) 
(4.23) 

(4.24a) 

(4.24b) 
(4.24c) 



The expressions for Cq and C\ were first obtained by Ferdinand and Fisher [[Tl| . The 
results for C 2 and C3 are new. In Table ^| we show the values of the coefficient C 2 {p) 
for several values of the aspect ratio p. 

Remarks 1. As shown by Ferdinand and Fisher using the Jacobi's transformations 
( |A.4| ), the coefficient Co(p) satisfies the identity 



C (p) = C (p 



-1^ 



TX 



logp 



(4.25) 



Indeed, from relations ( |4.21|) / ([4.23|) we conclude that these two coefficients have the 
right behavior under the transformation p —>■ 1/p [c.f.,( |3.25"a| )/(3.25b)1. Finally, we 
have tested numerically that 

CM = ^ (4.26) 

P 2 

is satisfied. This is a non-trivial test of the correctness of our result. 

2. We have also performed the following test: we have defined the function equal 
to the exact value of the specific heat Cn, c (n, p) minus the first three terms of the 



expansion ( 4.16 ). Then, we have fitted the result to several Ansatze. In particular, 
we show in Table [| the results for the Ansatz D 2 /n 2 + D 3 /n 3 + D 4 /n 4 + D 5 /n 5 . We 
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C 2 (p) 



1 


0.097119896855337 


2 


-0.326865280829340 


3 


-0.748187561687100 


4 


-0.877385104391125 


5 


-0.809168407448959 


6 


-0.682469414146146 


7 


-0.567445479079586 


8 


-0.483401539198706 


9 


-0.428249598449714 


10 


-0.394311952593824 


15 


-0.351150319692501 


20 


-0.349140134316672 


00 


-0.349065850398866 



Table 3: Values of the coefficient C^p) from (|4.22|) for several values of the torus 
aspect ratio p. 



have used the same data as for the energy fits in Section |^. The agreement between 
the values D 2 , D 3 and the exact values (|4.22|) / (|4.23|) /( ^23D is very good. It is also 



interesting to note that D 5 (p) « —\^2B 5 (p), where B 5 (p) is the coefficient obtained 
in a similar fit to the energy (see Table 



_p 1 2 3 

D 2 0.097119896855 -0.3268652808 -0.7481875617 
D 3 0.2922941107 0.2605022 0.2167250 
D 4 0.014792 -1.0635 -1.8668 
D 5 1.0326 059 045 

Table 4: Fits of the function C# jC (n, p) — (8/7r)logn — Co(p) — Ci(p)/n [cf. 
(P^)/ (|O0D / (|OlD 1 to the Ansatz D 2 n~ 2 + D 3 n~ 3 + D 4 n~ 4 + D 5 n~ 5 for several 



values of the torus aspect ratio p. 

3. By inspection from Table ^ the function C 2 (p) should have a zero at a non- 
trivial value of the aspect ratio p m i n between 1 and 2. We have evaluated numerically 
that value 

p min w 1.33544086 (4.27) 

By ( ^26|) , there is another zero of C 2 {p) at p m ^ n w 0.74881639. 



4. In the limit p — > oo the coefficients C{ tend to the following limits: 

Cq(oo) = - log + lE--r) (4.28a) 

7T V TV 4 / 





Ci(oo) 


= 












(4.28b) 




C 2 (oo) 


7T 












(4.28c) 




C 3 (oo) 


= 












(4.28d) 


Thus, 


only the coefficients associated to 


even 


powers 


of n 


~ l survive 


in 


this limit. 



5 Further remarks and conclusions 

The computation of the finite-size corrections to the energy and specific heat shows 
that, through order n~ 3 : 

(a) All the corrections are integer powers of the quantity n~ l . In particular there 
are no multiplicative or additive logarithmic terms (except for the leading term 
in the specific heat). 

(b) In the energy density we only find odd powers of rT 1 . In particular, the correc- 
tions of order rT 2 and n~ 4 are absent in Q3.21| ). 



In the specific heat we find both even and odd powers of n~ l . But the co- 
efficients of the corrections corresponding to odd powers of n~ l in (|4.19|) are 
proportional to the corresponding coefficients in the energy expansion ( |3 . 2 1| ) : 
we found that d(p) = -y/2Ei(p) for i = 1, 3 [cf. flOIft /flOJPI, and the numer- 



ical test performed at the end of Sections || and f| shows that the same ratio 
holds for the next coefficients C$/E 5 « — \[2. 

The natural question is whether these observations are general features that hold to 
all orders in n _1 . In this section, we will try to answer those questions. 

Let us first analyze what happens to the energy. We first note that in the expansion 
( |3.8| ) only even powers of n~ l occur. This expansion can be done to any finite order we 
want. The errors coming from neglecting the second sum in ( |3.5|) and from extending 
the sums X)r=i ^° ^2T=i are a ^ mos t of order 0(n~^ M ^), and they can be made as 
small as we want by making M in the definition of s(n) ( |3.6| ) as large as we need. This 
means that in the series expansions of the quantities logP, (i — 1, . . . , 4) (|3.13| ) only 
even powers of n^ 1 appear. Secondly, we need to check that only even powers of n^ 1 
occur in the expansion of logPo ( |3.15| ). The argument is simple: all derivatives f^ of 
the function f(p) = log(sinp+A/l + sin 2 p) are integrable over the interval [0, it]. This 
implies that the Euler-MacLaurin formula ( |B.4j ) can be applied to any arbitrary order. 
As /(0) = f(ir), the correction of order n _1 vanishes and only corrections with even 
powers of n~ l can occur. The same conclusion applies to the expansion ( 3.15Q and to 



the ratios Ri ( ^.20[ ). Thus, from formula ( pj.3[ ) we immediately conclude that only odd 



powers of n^ 1 appear in the finite-size corrections to the critical energy density. This 
result generalizes point (b) above. In particular, no logarithmic corrections occur in 
this expansion at any order. 



19 



Remark. The authors of Ref. |^2| made an argument to explain the absence of the 



L 2 correction in the energy. They started with the Fortuin-Kasteleyn representation 
of the (/-state Potts model on a graph G (q = 2 corresponds to the Ising model) |29], [30 



Z G (K) = J2 v N ^ G 'k Nc{G,) , (5.1) 



G'CG 

where N b and N c are respectively the number of bonds and connected clusters of the 
spanning subgraph G' C G, and v = e K — 1. Then, they introduced the standard 
finite-size-scaling Ansatz for the free energy assuming there are no irrelevant operators 

& 

fa{K) = U(K) + -}—W{L]l v {K - K c )) . (5.2) 

Ll X Lly 

Here G is a square lattice of size L x x L y , f TCg (K) is the regular part of the free 
energy, v is the usual critical exponent, and W(x) is an analytic function at x = 0. 
They expanded W(x) around x = and computed the mean values of N c and N b at 
criticality (K = K c ). They found that 

(N c ) = n c L x L y + AL\ jv + B (5.3a) 

(N b ) = n b L x L y + CL]l v (5.3b) 



with no constant correction to (Nb) ■ The internal energy ( |2.7| ) is related linearly with 
N b ; thus, in the Ising model at criticality, we have E c = A + Ai(p)L~ 1 with no higher- 
order corrections in L~ x . They concluded that the L~ 2 correction to the critical energy 
should vanish. Indeed, this argument also implies the stronger result that there are 
no corrections of any kind beyond order L^ 1 - - a conclusion that is unfortunately 
false! So it is not clear that the Ansatz ( |5.2j ) (even if we include irrelevant operators) 
is enough to reproduce the right finite-size expansion of the internal energy of the 
critical two-dimensional Ising model ( |3.21| ). A detailed discussion of this point will 
be published elsewhere [^T| (See also |14|, Section 3] for preliminary results). 



The analysis of the specific heat is a little more involved. Using the same procedure 
as for the energy, we conclude that the fourth (~ R3/R) and the last (~ R\j ' R 2 ) terms 
of ( |4.1| ) provide corrections with even powers of n _1 . Furthermore, the sixth term of 
Tj) [~ R±/{Rn)\ will give only odd powers of n~ l . 

Let us see nor what happens to the terms Qij ( 4.2a| -d). The argument for Qi,4 



( |4.2d|) is quite simple: the expansion of this quantity as a power series of rir/n will 
only contain even powers of n~ l . Indeed, the errors coming from dropping the sum 
Ylr=s{n) ^ n ^ e beginning of the computation, and for extending the sum Ylr=i ^° 
X^li & t the end of the computation are both at most of order 0(x^ ( - M ^ 1 ') with the 
choice ( |3.6| ) for s(n). Thus, the fifth term in ( |4. 1| ) [~ ^2RiQi,i/R] will have only 
corrections with even powers of rC l . 

Finally, let us analyze the behavior of Qi,+ ( [4.21] ) . In the evaluation of Q\ + we 



had to apply the Euler-MacLaurin formula (|B.4j ) to the function f(p) = sin 2fc x p with 
k > 1. The derivatives of this function are always integrable over the interval [0, tt], 
thus the expansion ( |4.5|) only contains even powers of n~ l . In the evaluation of Q\ + 
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we apply the Euler-MacLaurin formula (|B.4| ) to the function f(p) = 1/sinp — 1/p + 
l/(p — 7r). This function and all its derivatives are integrable over the interval [0, n]. 
Thus, Eq. ( f4.6| ) can be generalized to 

{77. — 1 OO 1 

p=i ^ fe=i J 

This expansion contains in principle both even and odd powers of n" 1 . We now plug 
in the well-known result |28| 



±1 logL + , E+ l. ±^ (5 . 5 ) 

p=l ^ fc=l 

with L = n — 1, and expand the resulting factors log(l — 1/n), (1 — 1/n) -1 , and 
(1 — l/ri)~ 2fc in powers of n _1 . If we express the result as Y^ fc a^n~ k : , we immediately 
see that cui = 0. The expression for the coefficients of the odd powers of n~ l is given 
by 

m=l x ' 



To prove c^+i = 0, we can apply the general Euler-MacLaurin formula ( |B.1| ) to the 
function /(x) = x 2k with n = and m — 1: 

x dx h> -: rr^-; rr = — «2Ai+i (5.7) 

2 ^ (2m)!(2fc-2m + l)! + V ; 

m=l v ' v ' 

This implies that in the finite-size-scaling expansion of Qi,+ only even powers of rT 1 
occur. The same holds for Q\ (|4.2e|) . 



In summary, we have seen that in the finite-size-scaling expansion of the specific 
heat at criticality only integer powers of n" 1 appear [except of course for the leading 
term (8/tt) log n] . Furthermore, all contributions to the specific heat ( |4.1| ) give even 
powers of n~ l , except for one: namely, the sixth term in ( |4.1|) , which has the form 

2vg lim ^(r,n,p) coth 
n r-*o R(0,n,p) 

If we compare it to the expression for the energy ( |3.3j ) we conclude that the coeffi- 
cients associated to odd powers of n^ 1 in the energy and specific-heat expansions are 
proportional. In particular [c.f., ( |2.12|) / (|2.13| )1, and recalling that E 2 k = 0, we obtain 



E k (p) _ J -1/V2 for k odd . Q) 



Ck(p) I for A; even 



This generalizes the results (|4.21|) / (|4.23|) . Note in particular that the ratio E k (p)/Ck(p) 



is independent of the aspect ratio p. 
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Let us conclude with a brief discussion about the generality of the above results. It 
is clear that the leading-term coefficients E and C 00 in (|2.12| )/ (|2.13|) do not depend 



on the boundary conditions, as they are bulk quantities. However, the finite-size- 
scaling coefficients E k (p) and Cfc(p) with k > 1 are expected to depend in general on 
the boundary conditions of the system. In particular, Lu and Wu []3T| have obtained 
the finite-size expansion of the free energy for a critical Ising model with two different 
boundary conditions (namely, a Mobius strip and a Klein bottle). They found an 
expansion of the form 

ft \ f t fl\P) , h(P) , /r 1n \ 

fc{n, p) = fo + + — 5- + • • ■ (5.10) 

n n z 

where the coefficients /&(/?) with k = 1,2 do depend explicitly on the boundary 
conditions. 

On the other hand, the fact that the ratio Ek(p)/Ck(p) is a p- independent number 
suggests that it might be universal, i.e., independent of the details of the Hamiltonian. 
It would be very interesting to investigate this possibility. A first step in this direction 
has been achieved by Izmailian and Hu |32] who computed the finite-size expansion of 
the free energy per spin fjy and the inverse correlation length £^- for a critical Ising 
model on several N x oo lattices with periodic boundary conditions (namely, square, 
hexagonal and triangular). On each lattice they found expansions of the form 

oo „ 



fc=l 

oo 



S 1 Y,wh (B-Ub) 



fc=l 



where the coefficients b k and /& do depend on the lattice. However, their ratio bk/fk 
is the same for all three lattices and equal to bk/fk — (2 2fc — l)/(2 2fe_1 — 1). They 
also computed the corresponding expansions for a quantum spin chain belonging to 
the two-dimensional Ising model universality class. They found that the ratio bk/fk 
takes the same value as in the Ising case. This result supports the conjecture that 
this ratio is a universal quantity. 

A Theta functions 

In this appendix we gather all the definitions and properties of the Jacobi's 9- 
functions needed in this paper. We follow the notation of Ref. [P|, which was adapted 



from Whittaker and Watson |33] . We define the Jacobi ^-functions 0j at z — in the 
following way:Q 



9 2 = 9 2 {0, e- wp ) = 29 {) e- np l A \\ (l + e - 2r7Tp ) 2 (A. la) 



L In this particular case, #i(0, e 7rp ) = 0. 

22 



) 



r=l 
oo 

4 = 4 (O, e~ wp ) = O J] (1 - e- (2r - J -" 

r=l 

The function 9o(p) is denned as 

#o = e { P ) = n (1 - e- 2 ^) 

r=l 

and it satisfies the following identity 

#0 = e^ /12 [^020304 



1/3 



These 0-functions ( |A.1|) satisfy the Jacobi's imaginary transformation 

2 (O,e-^) = p 1 ^ (0, e"*') 
3 (O,e-^) = p 1/2 3 (O, e"*') 
4 (O,e-^) = p 1/2 ^ 2 (0, e ^) 



(A.lb) 
(A.lc) 



(A.2) 



(A.3) 



(A.4a) 
(A.4b) 
(A.4c) 



B The Euler-MacLaurin formula 

The Euler-MacLaurin formula (see e.g. |34], Appendix B]) is the main tool we need 
to compute asymptotic series. The general form of this formula is given by 



m— 1 



£/(*) 



k=n 



dxf(x) -~[f{m) - f{n)] 



N 



B, 



+ £|^[/ (2p ~ 1 V)-/ (2p - 1) (n)] 



+ 



(2iV + l)!7„ 



dxf (2N+1 \x)B 2N+1 (x-[x\) (B.l) 



where B n are the Bernoulli numbers and B n (x) are the Bernoulli polynomials defined 
by 



Bn(x) = J2 ( I ) B ^ H ~ k 
fc=0 ^ ' 



(B.2) 



In this paper we are mainly interested in sums of the form 

1 aL-l 
n=0 



(B.3) 
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where p = 2irn/L. The asymptotic expansion of the sum ( |B.3|) in the limit L — > oo 
with a fixed can be obtained from (|B.1| ): 



aL-l r 2iTa 



1 r dn I 

n=0 JU 



TV / v 2k 



x 2JV+1 ^2™ 
Z7T 



where 



(2JV + 1)! \L J J 2it 

Lp 



^ 2N+1 \p)B 2N+1 (p) (B.4) 



S n (p) = 5n ( ^ 

Z7T 



2^ 



(B.5) 



The expression ( |B.4| ) gives the asymptotic expansion of the sum ( [B.3| ) in powers of L x 
up to order L~ 2N if the last integral in (|B.4|) is finite (i.e., if f^ 2N+1 \p) is integrable 



in the interval [0, 2wr]). If /(0) = /(27ra), then only even powers of L x occur in the 
expansion 
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